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SYMBOLS 


area of portion i of interface 

speed of sound 

diameter of cylinder 

specific internal energy 

interface between impacting materials 

portion of interface inside boundary defined by points of intersection 
of cylinder's radial rarefaction with interface 

impulse given to target material in time t£ by pressures acting at 
portion i of interface 

constant 

length of cylinder 
mass of cylinder 

penetration (maximum depth of crater) measured from undisturbed 
target surface 

pressure 

rarefaction (head of rarefaction wave) 
radial coordinate 
shock wave 

time measured from instant of impact 

time for shock wave to reach rear face of impacting cylinder 
time at which axial rarefaction reaches interface 

velocity of one-dimensional shock wave relative to undisturbed medium 

mass velocity behind one-dimensional shock wave relative to 
undisturbed medium 

velocity of cylinder at impact 


axial coordinate 



0 dimensionless parameter defined by equation (9) 

p = n - 1 
p density 

Subscripts 

H Hugoniot state 

o initial state 

p projectile or cylinder material 

r radial direction 

t target material 

x axial direction 

y impact case of low-density cylinder 

a impact case of simulating projectile 
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SIMULATION OF METEOROID-VELOCITY IMPACT BY USE OF 


DENSE PROJECTILES 
By Robert H. Morrison 
Ames Research Center 


SUMMARY 


A technique is described for simulating the impacts of semi-infinite 
targets by low-density, low-fineness-ratio cylinders at meteoroid velocities 
by impacting high-density, low-velocity projectiles. Conditions for simula- 
tion are stated and embodied in a model forming the basis of the technique. 

The feasibility of the technique is experimentally investigated by simulating 
the impacts of semi-infinite aluminum targets by polyethylene plastic cylin- 
ders of various low-fineness ratios and the same velocity of 11.3 km/sec. 

The penetrations were simulated at lower velocities to within 10 percent by 
impacting aluminum, steel, nickel, and copper projectiles, but only to within 
30 percent by impacting platinum projectiles. In application of the tech- 
nique, the impacts of semi- infinite aluminum targets by cylinders of poly- 
ethylene and porous aluminum (0.44 g/cm 3 ) at velocities of 15.2 and 22 km/sec, 
respectively, were simulated. 


INTRODUCTION 


The simulation of meteoroid impact has, for the most part, been beyond 
the capabilities of light-gas-gun facilities. The highest impact velocity 
attained in these facilities, 11.3 km/sec, is just at the lower end of the 
meteoroid velocity spectrum of 11 to 73 km/sec (ref. 1) . 

This velocity (ref. 2) was obtained with a low-density, plastic cylinder 
that impacted an aluminum target. Calculations indicate that the initial 
pressure produced in the target material by this impact is less than that 
produced by much slower, higher density projectiles. This observation sug- 
gests that with the proper dimensions and velocities, higher density projec- 
tiles could generate pressure pulses in the target material which would 
approximate that generated by the plastic cylinder impacting at a velocity 
exceeding 11.3 km/sec, and the projectiles thereby could simulate, to some 
degree, its impact damage at presently unattainable velocities. In the same 
manner, the impacts of low-density meteoroids at higher velocities could be 
simulated. 

As a step toward this goal, design criteria are herein proposed for the 
dimensions and velocities of dense projectiles to simulate the impacts of 
low-density, low-fineness-ratio cylinders at meteoroid velocities. The model 
in which these criteria are given and the results of experiments designed to 



test the technique are presented... Penetration data are included for effec- 
tively semi- infinite 2024-T351 aluminum targets impacted by aluminum, steel, 
nickel, copper, and platinum cylinders at velocities ranging from 4.41 to 
8.35 km/sec. 


THEORETICAL CONSIDERATIONS 


Early Stages of Hypervelocity Impact 

During the early stages of hypervelocity impact, much of the process 
can be analyzed in terms of highly simplified flows: one-<Jimensiqnal flow 

behind a normal shock and expansion-wave heads propagating from free surfaces 
through material at uniform pressure. If enough of the impulse given to the, 
target occurs during this '.'simple" .part? of the flow, it should be possible 
to define conditions for impact equivalence between projectiles of widely 
different physical properties. 

, The axisymmetric, hypervelqcity impact of a semi-infinite, target by ,a. low- 
fineness-ratio cylinder is depicted at the instant of impact in. figure 1. The 
cylinder's velocity v, diameter d, length t, and initial density p 0 p, >: 
which in the general case is different, from .the initial density of the target, 
p ot , are indicated. The orthogonal coordinate system x - r is fixed in the 
target with the x axis coinciding with the axis of the cylinder. From the 
instant of impact, the materials of both the cylinder and the target willbe 
compressed to very high pressures by shock waves. These pressures, greatly 
exceeding the strengths of both materials, will start a nonsteady, three- 
dimensional compressible flow. 

• The salient features of this compressible flow at two early times after 
impact ate, shown ; in figure 2, for a . ease where . the speed of sound at the 
initial impact (Hugoniot) pressure in the cylinder material is less than that 
of the target material. , 


In figure 2 (a) the wave pattern of this flow in a cross-sectional plane 
through the x. axis has been estimated for a time soon after impact., The 
shock waves Sp and S t . are shown propagating away from the interface I 
between the materials of the cylinder and target, compressing these respec- . 
tive materials. Also shown are the heads R r p and R r t of the radial rare- 
faction waves, which start inward at the instant of impact from the surface 
of the cylinder. Since the speed of sound in the target material is greater, 
the radial rarefaction in this material precedes that in the cylinder material 
along the interface. Because pressure waves are transmitted across the . 
interface, expansion Mach waves are formed in the cylinder material. However, 
in figure 2, the cylinder's radial rarefaction is shown as though its propa- 
gation were unaffected by these expansions. When the speed of sound at 
impact is greater in the cylinder material, Mach waves are generated in the 
target material instead. (In similar materials, no such waves are formed.) 

A radial rarefaction and the associated Mach waves form a boundary, inside of 
which the flow, is still one- dimensional and the two shock waves and the 
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interface are planar. - Outside this boundary,, the flow is- three-dimensional 
and the shock waves and interface are curved. < • / •. 


At a later time, the shock Sp will; have reflected from the cylinder's 
reap face as a rarefaction'. The wave .pattern at this time is estimated in 
figure. 2(b) . Shown again- are 'the radial rarefactions in the target material, 
which by this time are , approaching the ; x- ' axis , . having .penetrated most of 
the one- dimensional flow,. Inside the region bounded by the radial expansion, 
system,, the axial rarefaction head. R X p remains planar, 

... The material and wave motions occurring along the- x axis during these 
early stages of the impact are better, shown- in the waye diagram .. of figure, 3.. 
The shock wave Sp propagates through the cylinder and, upon encountering 
the cylinder's rear face at time tj, , is reflected as the axial rarefaction 
R X p, which starts the isentropic release of the material from the Hugoniot 
pressure. Subsequently, at time t 2 w the axial rarefaction overtakes the 
interface I. For the fineness ratios being considered, the axial rarefac- 
tion arrives at this point before either of the radial rarefactions . 

• The cylinder velocity v and the, wave - arrival ..times tp and t 2 . may 
all be related to the impacting, cylinder.'s dimensions- and the one-dimensional 
flow parameters (see appendix A) . These parameters are the speeds U and u, 
relative' to the undisturbed medium,; of the shock waye and- the . compressed 
material,, respectively, and ..the .speed of sound c H - an this compressed: mate.-, 
rial. at pressure p^.. These parameters- may he calculated, for 1 various, mate.-’ 
rials and any given pressure pjj,,as indicated in appendix B. The relations 
are . • _ > 

v = Up + u t (1) 

t! = l/V v . (2) 

t 2 = Kl/U^ l/n Hp c Hp ) (3) 

where 

n R = U/CU, - u) 

u, U, and c^ must be evaluated for the materials of the 
p) and the target (subscript t) at the same Hugoniot 


Simulation Model * " ' ' 

It is postulated that the initial part of. the pressure pulse generated 
in a semi-infinite target by the hypervelocity' impact of a low-density, low- 
fineness-ratio cylinder can be-, simulated and that this impact can thereby be 
simulated with the impact of a slower, higher density cylindrical projectile 
if the following conditions are met. The Hugoniot pressure p^ and the 


The flow parameters 
cylinder (subscript 
pressure. 
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arrival time t 2 of the axial rarefaction at the interface are to be the 
same for both cases. Furthermore, the same impulse J is to be given to 
the target material in this time by pressures acting at the portion i of 
the interface inside the boundary that is delineated by the points of inter- 
section of the cylinder's radial rarefaction with this interface. To delin- 
eate this boundary, the assumptions are made that the precedence of a target's 
radial rarefaction causes a negligible reduction in the pressure at the inter- 
face from the Hugoniot pressure ppj before the arrival of the cylinder's 
radial rarefaction, and that this precedence therefore has a negligible effect 
on the latter wave's propagation. Thus, the portion i of the interface is 
approximately planar and at a constant pressure p H . The above conditions for 
simulation may be written, respectively, as 

PHa = ?Hy 
*2 a ~ *2y 

J = J 
o y 

where the subscripts y and a denote the impact cases of the lower density 
cylinder and the simulating projectile, respectively. 


(4) 

(5 a) 
(6a) 


The preceding conditions and assumptions are sufficient to specify the 
velocity, length, and fineness ratio of a simulating projectile in terms of 
these same quantities for the low-density cylinder and the one-dimensional 
flow parameters. The velocity of this projectile is derived by substituting 
equation (1) into the equation 


u 


to 


ty 


(7a) 


which is implied by the condition of equation (4) . ■ The velocity thus 
obtained is 


v = v + u - u 
a y per py 


(7b) 


The length of this projectile is obtained from equations (5a) and (3) as 


l 

a 


K - »'VH»p) u . 

( 1/U p * lv 


(5b) 


Its fineness ratio is derived from equation (6a) by first developing an 
expression for the impulse J as follows. If is the area of the portion 

i of the interface, the impulse J may be expressed as 


J = 



dt 



dt 


However, the area is given by 



A i - 7 ( d - 2c Hp *)* 


( 8 ) 


Substituting this expression for into the one for J, performing the 

integration, and rearranging the resultant expression yields 

J = 7 ^H t 2 c Hp® 

where 0 is a dimensionless parameter defined as 


0 


1 

_ c Hp( 1/U p * 1/ "Hp c Hp)( J/d ) 



(9) 


If the latter expression for J is substituted and the conditions of equa- 
tions (4) and (5a) are applied, equation (6a) becomes 


~ ( c Hpy^ c Hpcr) 9 y 


(6b) 


By use of the identity of equation (9), the projectile's fineness ratio is 
found to be 

- [ C Hp 0 ( 1/U P +1 ^' 1 Hp C Hp) a (' * v /5 r r T7l)j 1 


( 10 ) 


where 0 is related to the low-density cylinder's properties by 
equation 0 (6b) . 

In summary, the velocity, length, and fineness ratio of a cylindrical 
projectile to simulate the impact of a low-density, low-fineness -ratio 
cylinder are given, respectively, by the equations 
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and 


h H = u/(u - u) 


■ EXPERIMENTS 


The experimental investigation .consisted of two parts... First, the 
applicability of the preceding simulation relationships was determined by 
simulating the impact of the previously mentioned plastic cylinder of ref- 
erence 2. The cylinder wais polyethylene with a density jof 0.95 g/cm 3 and a 
fineness ratio of 1/3. " To examine somewhat -more generally- the applicability 
of the simulation technique, the impacts of. polyethylene cylinders with lower 
fineness ratios and the same velocity (11.3 km/s ! ec) were also simulated. In 
the second part of the experimental investigation, an impact of the 1/3- 
fineness -ratio polyethylene cylinder Was simulated at a higher: velocity of' 
15.2 km/sec. For each simulatioh,'. the 'penetration (target-crater depth) so 
obtained was compared with that estimated from 

P/d = 0.34(//d) l/4 v 2/3 (11) 

(ref. 3), which is applicable to the hyperveiocity impact of 2024-T351 alumi- 
num by polyethylene cylinders with fineness ratios, from 1/6 to 1 . 

To simulate the impacts in the first part of the experimental investiga- 
tion, the following projectile materials were used: 2024-T351 aluminum, AISI 

1018 steel, pure nickel , OFHC ‘ copper , and pure platinum.' For the second part, 
only OFHC copper was used. The velocities of the various simulating- projec- 
tiles for yielding Hugoniot pressures of 1.05 and 1.70 megabars and thus 
simulating/polyethylene impacting at 11.3 and 15 j 1 -km/See, respectively; 'are 
listed in' -table 1. The velociti'es were calculated for eadh pressure from 
equation (7b) by substituting the -velocity v of polyethylene for v and the 

values calculated as indicated in appendix B for u and u . fJ 

■ pa PM 

Projectiles of these materials were machined with diameters near 5 mm 
and various low-fineness ratios. The diameters and lengths were measured to 
within 0.06 and 0.6 percent, respectively. The projectiles were then 
launched, sabot-mounted, at the required velocities- in the ballistic ranges 
employed in the investigations of references 2 . and 4. The light gas guns of 
these ranges were equipped with rifle barrels, which imparted spin to the 
projectiles to ensure attitude near zero angle of attack, and separation of 
the segmented sabOjts. . Spark shadowgraphs taken at -.numerous stations on these 
ranges provided orthogonal views^' of the- projectiles in flight and thereby 
gave condition and attitude data. The angles of attack did not exceed 5° 
and therefore were in the range found in reference 3, for which any anomalies 
due to angles: of attack are outweighed by the inherent variation in crater 
formation. These shadowgraphs, together with the measurements, for intervals 
of time and distance : between the stations y made possible, in each case a 
determination! of the impact velocity to within 1 percept. ■' 
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The targets were of the same material (2024-T351 aluminum) and size 
(10 x 10 x 5.1 cm) as those of references 2 and 3'. The depths of the resul- 
tant craters were measured relative to the undisturbed target face to within 
0.6 percent by an especially adapted dial indicator. 

RESULTS AND' DISCUSSION - ‘ \ 

The pertinent data and the results of the simulation tests are given in 
tables 2 and 3, respectively. In table 3, the dimensionless penetration 
(P/d) corresponding to the velocity v a , together with the density P ocr , 
length l a , and fineness ratio (Z/d) 0 , is listed for each simulating projec- 
tile. The velocity v y , length 7 y , and fineness ratio (Z/ d) y of the poly- 
ethylene cylinder whose penetration P y was Simulated by each projectile and 
the penetration of this projectile normalized with respect to that estimated 
for the corresponding polyethylene cylinder P 0 /P y - are also listed, . ’In addi- 
tion, the maximum deviation of tbe dimensionless penetration P/d is given 
for each projectile from a curve described by' an equation of tbe form 
P/d = Kv 2 / 3 and fittpd to the data by the method of least squares. This 
relation between penetration and impact velocity was shown to apply for 
hypervelocity impact tip to the 11.3 km/sec in reference 1 2, The maximum' 
deviation of P/d, within 3 percent for each case where several data points 
were obtained, shows that this relation adequately describes the data over' the 
velocity range tested. In each case,, the equation was used to estimate (P/d) a 
at the simulating velocity v 0 . The length l y and (Z/d) y were calculated 
by means of equations (5b) and (10), respectively j' where'aS the normdlized' 
penetration was calculated from the' equation ' 

P„ - 1 (l/’Jp * (i’/di, 

P„ "OT-U/Up ;■» ■ (t/,d) d; Vy2/ 3- ■■ 

which is derived by substituting equations. (5b) and (11) into the identity . 

p„ (p/d). c/iip v \ - ^ ' ■ 

....... : p M; ( p /d) u .cz/dio z u . .; :..j 

In using equation (12) for the smallest fineness ratio aluminum projectile 
and its corresponding O’. 118-fineness ' ratio polyethylene cylinder, it was • 
assumed that the l/d . relation of equation (11) Still- held-, ' 

For each simulating projectile, the normalized penetration P 0 /P y is 
plotted Versus the fineness ratio (Z/d) y of the corr espondihg ; polyethylene ' 
cylinder in figure 4. This figure and table 3 show that the aluminum, s te'el) 
nickel, and copper projectiles, ' having fiheness ratios 'from* 0 ,0816' to 0, 224 > 
simulated with good accuracy (to within 10 percent) the penetrations of their 
corresponding polyethylene cylinders, having fineness ratios from 0.118 to 
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0.313. However, poor results were obtained with the platinum projectiles, 
since the penetrations were deeper than those of their corresponding 
polyethylene cylinders by as much as 30 percent. 

The impacts of the steel, nickel, copper, and platinum simulating pro- 
jectiles correspond to the previously discussed case of figure 2, where the 
speed of sound in the target material at high pressure is greater than that 
of the cylinder material, and therefore the radial rarefaction in the target 
precedes that in the cylinder along the interface. For the impacts of the 
corresponding polyethylene cylinders, the opposite is true. The fineness 
ratio of each simulating projectile and its corresponding polyethylene cyl- 
inder, however, was sufficiently low that the axial rarefaction reached the 
point on the axis at the interface before either of the radial rarefactions . 

The poor results obtained with the platinum projectiles cannot be attrib- 
uted solely to experimental errors since these errors were no larger than 
those for the smallest fineness ratio aluminum projectile, which simulated 
closely the polyethylene impact. Nor is the variation in crater depths, as 
represented by the maximum deviation of P/d from the fitted curve, large 
enough to account for the 20 to 30 percent deeper penetration of platinum 
since this deviation was only 1 percent. Although a different equation of 
state was used in calculating the one-dimensional flow parameters for plati- 
num (see appendix B) , its use also is not thought to be solely responsible 
for the much deeper penetrations . 

Rather, the poor results obtained with the platinum projectiles are 
probably due to the radial and axial expansions of projectile material being 
much slower in each of these cases subsequent to time t 2 - The slower radial 
expansion is evident from table 4 where it is shown that, despite selection 
of fineness ratio and length to duplicate the impulse J of the impact of 
the 1/3-fineness-ratio polyethylene cylinder, the area of the one-dimensional 
flow at time t£ is much larger than in the polyethylene case. The same is 
true for simulation with copper and aluminum, but to a lesser extent. The 
slower axial expansion, on the other hand, is evident if one considers the 
reflection of the axial rarefaction from the interface, which must occur for 
impabts of dissimilar materials because of the acoustic mismatch present 
there. The much greater acoustic impedance of the platinum causes a large 
negative reflection of the axial rarefaction, which, in turn, causes the ini- 
tial pressure decay at the portion i of the interface to be much slower than 
for the polyethylene case where a positive reflection occurs. The pressure 
decay is also slower for simulations with copper and aluminum, but again to a 
lesser extent. The much slower radial and axial expansions in the platinum 
cause the shock pressure in the aluminum target to be sustained to greater 
depths, thus producing greater penetrations. This suggests the inadequacy of 
the simple criteria for projectile dimensions when the physical properties of 
the projectile materials differ vastly from those of the low-density cylinder. 
However, the use of such projectile materials in simulating actual meteoroid 
impacts to develop a design criterion would result in an overestimation of 
meteoroid penetration ability and, therefore, lead to a conservative design 
for meteoroid protection. 
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Profiles of the craters produced by the simulating projectiles impacting 
at Hugoniot pressures near 1.05 megabar are depicted in figure 5. As can be 
seen from the photographs, the crater shapes are more conical than hemispher- 
ical, which was the result obtained in reference 2 for the impact of the high 
speed polyethylene cylinder (see fig. 6). However, in comparison, the cra- 
ters made by the platinum projectiles have relatively flat bottoms. The 
craters for the aluminum, copper, and platinum projectiles of fineness ratios 
0.224, 0.151, and 0.122, respectively, (see figs. 5(a), (d) , and (e)) nearly 
simulate the impact of the 1/3-fineness-ratio polyethylene cylinder of ref- 
erence 2. Their profiles are compared in figure 7. Except for the platinum 
impact crater, these craters are very similar to that produced by the poly- 
ethylene cylinder. 

The results of the tests indicate the feasibility of using slower moving 
but more dense projectiles to simulate the impacts of semi-infinite targets 
by low-density, low-fineness-ratio cylinders at meteoroid velocities. It is 
emphasized that the technique is restricted to impacts on semi-infinite tar- 
gets, but an analogous treatment of thin-target impact might prove useful. 

In particular, the technique can be used for simulating the impacts of low- 
fineness-ratio cylinders of polyethylene and porous aluminum at presently 
unattainable meteoroid velocities. 

In this application of the technique, the impact of a 1/3-fineness-ratio 
polyethylene cylinder at 15.2 km/sec was simulated by a copper projectile. 

The results are given in table 3 and figure 8. The simulated penetration 
agrees with that estimated by equation (11) to the same accuracy as that 
obtained for the lower velocity simulation. This agreement implies that the 
velocity exponent of this equation applies to impacts of polyethylene at 
velocities up to 15 km/sec. It is interesting that the shape of the simulated 
impact crater in figure 8 is more hemispherical than conical, which is a 
result just opposite to that obtained for the lower velocity simulation. 

As a further application of the technique, the impacts of porous alumi- 
num cylinders with density of 0.44 g/cm 3 that are simulated by (1) the high- 
velocity polyethylene cylinder of reference 2 and (2) the highest velocity 
copper projectile of table 3 were determined. The impacts of particles with 
this density are of interest since various astronomers have suggested densi- 
ties of about this value for cometary meteoroids (ref. 5) which are thought to 
be porous. The results are given in table 5, which lists the velocity Vy, 
length Z.y, fineness ratio (7/ d)y, the dimensionless penetration (P/d)y, 
and the ratio of the penetration to the cube root of the mass (P/m 1 / 3 ^ for 
the simulated porous aluminum cylinders. These quantities were calculated 
from the data for the simulating projectiles by using the values of the one- 
dimensional flow parameters for the porous aluminum as obtained from the 
appendixes of references 6 and 7. Also listed in table 5 for comparison is 
the ratio (P/m 1 ' 3 )^ for the impact of semi-infinite aluminum by an aluminum 
cylinder with the same porosity and velocity as the simulated cylinders, but 
with a fineness ratio of unity. The ratio in this case was determined from 
the theoretical results of reference 6. 
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The velocities of the porous aluminum projectiles, the impacts of which 
were simulated by the referenced polyethylene cylinder and the copper projec- 
tile, represent gains of 59 and 95 percent, respectively, over the velocity 
of the referenced cylinder. However, the penetration simulated by copper 
probably overestimates that which would actually be produced by the porous 
aluminum cylinder since the physical properties of these two materials differ 
vastly. In fact, the penetration simulated by copper is about 10 percent 
greater than that predicted by extrapolating from the penetration simulated 
by polyethylene according to the relation P/d « v 2 ' 3 . Nevertheless, using 
these simulated penetrations in determining required meteoroid protection for 
spacecraft should lead to a conservative design. 

The ratios (P/m 1 / 3 )y for the simulated porous aluminum cylinders do 
not necessarily apply to impacts of cylinders of other fineness ratios, but 
these ratios do agree rather well with those obtained theoretically for the 
unit-fineness-ratio cylinders. The maximum difference is only about 10 per- 
cent, which was obtained for the impact simulated by copper. If, for the 
porous aluminum cylinders, (P/m 1 / 3 )^ is indeed not a function of fineness 
ratio, then there is no shape effect for these cylinders similar to that 
found for the impacts of polyethylene cylinders in reference 3. 


CONCLUSIONS 


The results of the tests indicate that the penetrations into semi- 
infinite aluminum targets by polyethylene cylinders of various low-fineness 
ratios impacting at 11.3 km/sec were simulated to within 10 percent by the 
impacts at lower velocities of aluminum, steel, nickel, and copper projec- 
tiles, but only 'to within 30 percent by the impacts of platinum projectiles. 
From these results it is concluded that the technique is feasible for simu- 
lating the impacts of semi-infinite targets by low-density, low-fineness- 
ratio cylinders at meteoroid velocities, within certain limits. The 
principal limitation seems to be that the sound speed in the simulating 
projectile cannot be too small compared to the sound speed in the cylinder, 
the impact of which is to be simulated. The impacts of semi-infinite 
2024-T351 aluminum targets by cylinders of polyethylene and porous aluminum 
(0.44 g/cm 3 ) at velocities of 15.2 km/sec and 22 km/sec, respectively, were 
simulated by this technique. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., 94035, Nov. 5, 1969 
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APPENDIX A 


DERIVATIONS 


The velocity v can be related to the one-dimensional flow parameters if 
the boundary conditions of equal pressure and velocity are applied on either 
side of the planar interface (see fig. 2) . These conditions are 


and 


fH P * Phi = % (A1) 

v - Up = u t (A2) 


The first condition requires that the mass velocities Up and ut be evaluated 
at the same pressure ppj, and the second yields the relation 

v = Up + u t (A3) 

for the velocity. 

The relations for the arrival times t^ and t 2 can be derived in a 
manner similar to that indicated in reference 8. 
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APPENDIX B 


CALCULATION OF ONE-DIMENSIONAL FLOW PARAMETERS 


The parameters of the one-dimensional flow behind a planar shock wave 
are obtained by application of the conservation laws of mass, momentum, and 
energy across this wave. By choosing a coordinate system fixed in the 
undisturbed medium, one may derive the familiar Rankine-Hugoniot equations 

n H = p h /p o = u/(u - u) (Bl) 


P H - Pq = p o uU ( B 2a) 


E H - E o ■ (ph * Po)( 1/O o - 1 /p hV 2 

where U and u are the speeds, relative to the undisturbed medium, of the 
shock wave and the shock-compressed material, respectively. The symbols 
p, p, and E represent, respectively, the pressure, density, and specific 
internal energy of the material at the initial state (subscript o) and at 
the shocked Hugoniot state (subscript H) . For convenience, E Q may be 
chosen to be zero; and since for hypervelocity impact p^ >> p Q , the last 
two equations may be approximated as 


PH * l’o uU 

(B2b) 

Ph "h - 1 

^Pr, D u 

(B3b) 


where the definition of equation (Bl) has been substituted. 


To solve for the four unknowns n^, Ph> u, and U, a fourth equation, 
the equation-of-state of the material, is required. In this report, two 
forms for this equation were used for nonporous materials. The first of 
these is the semiempirical equation of Tillotson (ref. 9), 


P (p ,E) = (a + 1 ) Ep + Ay + By 2 (B4) 

\ E/E 0 n + 1/ 

where n = p/p Q , y = q - 1, and a, b. A, B, and E Q are constants. The 
values of these constants for various materials are listed in table 6, as 
given by or adapted from references 9 and 10. This formulation represents 
a best-fit interpolation between calculations from Thomas -Fermi -Dirac theory 
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for very high pressures and experimental data at low pressures and is 
estimated to be accurate to within 5 percent of the Hugoniot pressure for 
pressures below 5 megabars. The second form of the equation of state (Enig, 
ref. 11) is based on the "p-u mirror- image" approximation. This equation, 

4Sp(p.u H ) = (4Sp H ♦ p o C2)exp[-4Sp o (l/p - 1/p h )] - CBS) 

(the notation for S in ref. 11 is a) makes use of the linear approximation 

U = C 0 + Su (B6) 

(C 0 and S are empirical constants) , which is applicable for many materials 
over a large pressure range. The values of these constants for various mate- 
rials and the range of pressures for which they are applicable, as obtained 
from reference 12, are listed in table 7. 

The flow parameters may all be expressed in terms of the Hugoniot pres- 
sure ph by equations, the convenience of which is determined by the 
equation of state used. If Tillotson's equation is used, the most convenient 
equations for these parameters are 


u 



(B7a) 


U = 



(B8a) 


C H = fl (PH' "h) 



3pgp H /4PoEo * 4 (l * P H /2) 

(Vh / 2 o o E o ' "h ) 2 - 


(B9a) 


where c^ is the speed of sound in the shocked medium at the Hugoniot pres- 
sure. The first two equations are obtained by solving equations (Bl) and 
(B2b) . The last equation is found by applying the definition of the speed of 
sound, c 2 = (8p/9p)g, to equation (B4) and substituting into the result 
equation (B3b) and the thermodynamic relation (3E/9p)s = p/p 2 . The Hugoniot 
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pressure p^ in these equations is given by equation (B4) as 
PH = PhCph>%) • However, if equations (Bl) and (B3b) are substituted into 
this equation for p H , riu is obtained implicitly as = ^(p^) • Therefore, 
the flow parameters are functions only of the Hugoniot pressure 


u = u(p H ) 

(B7b) 

U = U(p H ) 

(B8b) 

C H = c h(Ph) 

(B9b) 

If Enig's equation is used, the most convenient equations 
are 

for these parameters 

u = -C 0 /2S + / (C 0 /2S) 2 + p H /p 0 S 

(B10) 

U = Cq + Su 

(B6) 

c H = (U + Su) (U - u)/U 

(Bll) 


Equation (BIO) is obtained by solving equations (B2b) and (B6) , whereas equa- 
tion (Bll) is obtained by applying the definition of c to equation (B5) and 
substituting into the result equations (Bl) and (BIO) . 

For the simulation tests, the one-dimensional flow parameters were 
calculated by using Tillotson's equation for polyethylene, 2024-T351 aluminum, 
AISI 1018 steel, nickel, and OFHC copper and by using Enig's equation for 
platinum. For porous aluminum of density of 0.44 g/cm 3 , the values of these 
parameters were obtained from the appendixes of references 6 and 7. The use 
of Enig's equation for platinum was necessary since the constants of 
Tillotson's equation for this material were not available. The constants 
used for polyethylene were adapted from those given in reference 10 for a 
slightly lower density polyethylene arid are listed in table 6. 

In addition, it was assumed that the equations of state for aluminum 
and steel alloys could be approximated by those for pure aluminum and pure 
iron, respectively. This approximation is good in the case of the aluminum 
alloy, since the values of u, U, and c H calculated for this material at 
Pjj = 1.000 megabar by using Enig's equation are within 2 percent of the cor- 
responding values calculated by using both this equation and that of Tillotson 
for pure aluminum. 

Using Enig's equation to calculate the one- dimensional flow parameters 
for Al, Fe, Ni, Cu, and W at p^ = 1.000 megabar gives values that agree well 
with those calculated by using Tillotson's equation. For each of these mate- 
rials, the values of u and the values of c H , calculated by using the two 
different equations of state, agree to within 1 percent and 3 percent, 
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respectively. Furthermore, the values of U agree to within 1 percent for 
each of these materials except Fe, for which the values agree to within only 
6 percent. The agreement of these values provides confidence in the values 
of the parameters for Pt, for which only Enig's equation was used. 
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TABLE 1.- SIMULATING PROJECTILE VELOCITIES 
FOR IMPACTS ON ALUMINUM TARGETS 


Projectile 

material 

Velocity, v , 
km/sec, at - 

p u = 1. 05 Mb 

ri 

P H = 1.70 Mb 

Polyethylene 

=V y =11.3 

=V y =15 .2 

2024-T351 A1 

7.40 

10.3 

AISI 1018 steel 

5.64 

7.83 

Ni 

5.39 

7.53 

OFHC Cu 

5.48 

7.61 

Pt 

4.67 

6.50 

1 




TABLE 2.- SUMMARY OF PERTINENT EXPERIMENTAL DATA 


Round 

Projectile 

Penetration, 

P, 

mm 

Projectile 

Material 

Density, 

g /2£3 

Diameter, 

d, 

mm 

Length , 

i, 

mm 

Mass, 

m, 

g 

Velocity, 

v, 

km/sec 

Mean 
length , 
mm 

Mean 

fineness 

ratio 

568 * 

Polyethylene 

0.95 

5.702 

1 . 

8 

0.0453 

11.3 

7.54 

1.8 

0.32 

989 

2024-T351 A1 

2. 

77 

5.075 

.411 

.0233 

7.40 

5.23 

.4 

14 

.0816 

994 





5.075 

.417 

.0233 

5.96 

4.57 





99S 





5.072 

.417 

.0232 

8.19 

5.56 





997 





5.075 

.411 

.0230 

8.35 

5.69 





990 





5.072 

.818 

.0467 

7.55 

6.50 

.822 

.162 

991 





5.075 

.826 

.0469 

5.86 

5.44 

.822 

.162 

J-108 





5.080 

1 . 

10 

.0616 

6.56 

6.38 

1.10 

.216 

J-117 





5.090 

1 . 

09 

.0611 

6.76 

6.55 





J-119 





5.085 

1 . 

09 

.0605 

6.95 

6.81 





J-122 





5.080 

1.10 

.0615 

6.75 

6.68 





J-123 









.0613 

6.97 

6.73 





J-124 









.0616 

6.61 

6.40 





J-125 









.0610 

7.27 

6.96 





J-126 









.0610 

7.22 

6.83 





J-127 





. 




.0614 

6.85 

6.66 



1 


J-203 



1 


5.362 

1 - 

20 

.0752 

6.70 

6.88 

1.20 

.224 

J-114 

AISI 1018 Steel 

7. 

84 

5.080 

.686 

.1079 

5.05 

7.11 

.687 

.136 

J-115 


7. 

84 



.688 

.1075 

5.31 

7.29 

.687 



J-116 

Ni 

8.86 



.691 

.1209 

5.06 

7.49 

.691 

! 


J-197 

OFHC Cu 

8.90 



.460 

.0812 

5.17 

6.30 

.460 

.0905 

J-213 







.518 

.0917 

5.44 

6.83 

.526 

.104 

J-214 





J 


.533 

.0953 

5.56 

7.14 

.526 

.104 

J-109 





i 


.645 

.1122 

5.80 

7.87 

.654 

.129 

J-110 





5.105 

.650 

.1102 

6.03 

7.98 





J-lll 





5.080 

.655 

.1088 

5.85 

7.75 





J-112 







.665 

.1123 

5.28 

7.24 





J-259 





< 


.655 

.1178 

6.06 

8.28 



1 


J-198 





4.821 

. 

729 

.1180 

5.42 

7.75 

.735 

.152 

J-205 







. 

744 

.1208 

5.56 

7.92 





J-206 







. 

739 

.1201 

5.74 

8.08 





SR- 196 







.732 

.1190 

6.89 

9.27 





SR- 197 





4.811 

. 

732 

.1181 

7.20 

9.52 





SR- 198 





4.821 

. 

742 

.1200 

7.63 

9.88 





SR-200 





4.811 

.726 

.1171 

7.82 

9.88 





J-211 

p.t 

21.42 

5.090 

.483 

.2098 

4.72 

9.35 

.481 

.0945 

J-249 







.475 

.2048 

4.62 

9.02 





J-251 







.485 

.2086 

4.70 

9.27 





J-207 





4.559 

• 

551 

.1923 

5.37 

10.2 

.550 

.121 

J-209 







.546 

.1912 

4.41 

9.12 





J-210 







.554 

.1935 

4.69 

9.42 



\ 



*From reference 2 
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TABLE 3.- RESULTS OF SIMULATION TESTS 
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TABLE 5.- SIMULATED IMPACTS OF POROUS ALUMINUM (0.44 g/cm 3 ) CYLINDERS 


Simulating projectile 



Porous aluminum cylinder 



Velocity, 

Length , 

Fineness 

Velocity, 

Length , 

Fineness 

Dimensionless 

CP/>n 1/3 ) y , 

(P/m l/3 ) u for 

Material 

v 


ratio. 

v u . 

t 

ratio. 

penetration. 

(Z/d) = 1, 


km/sec 

mm 

(Z/d) 0 

km/sec 

mm 

a/ d)„ 

0Vd)y 

cm/g 1/ 3 

cm/g 1/3 

Polyeth- 

11.3 

1-8 

0.32 

18 

2.5 

0.47 

1.4 

2.6 

2 . 6 a 

ylene 
OFHC Cu 

7.61 

0 v 735 

.152 

22 

2.6 

.47 

1.8 

3.3 

3 . 0 a 


From theoretical results of reference 6 


TABLE 6.- CONSTANTS FOR TILLOTSON'S EQUATION OF STATE 


— 

Material 

a 

b 

A, 

Mb 

B, 

Mb 

E o» 

Mb-cm 3 /g 

Normal density, 
p 0 » g/ cn > 3 

Sources 

Polyethylene 

0.6 

2.0 

0.075 

0.02 

0.07 

0.92 

References 9 and 10 

W 

.5 

1.04 

3.08 

2.50 

.225 

19.17 



Cu 

.5 

1.5 

1.39 

1.10 

.325 

8.90 



Fe 

.5 

1.5 

1.279 

1.05 

.095 

7.86 



A1 

.5 

1.63 

.752 

.65 

.050 

2.70 



Be 

.55 

.62 

1.1734 

.55 

.175 

1.85 



Ti 

.5 

.60 

1.03 

.50 

.070 

4.51 



Ni. 

.5 

1.33 

1.912 

1.50 

.090 

8.86 



Mb 

.5 

1.02 

2.713 

1.65 

.045 

10.20 



Th 

.4 

.86 

.531 

.50 

.025 

11.68 



Polyethylene 

.6 

2.0 

.075 

.002 

.07 

.95 

Assumed 


TABLE 7.- EMPIRICAL CONSTANTS IN THE EQUATION U = C 0 + Su 

[From ref. 12] 


Material 

Normal density, 
p 0> g/cm 3 

C Q , km/sec 

S 

Applicable pressure range, 
Mb 

W 

19.17 

4.005 

1.268 

0.395 - 2.074 

Cu 

8.90 

3.958 

1.497 

.883 - 1.444 

Fe 

7.86 

3.768 

1.655 

.43 - 1,967 

A1 

2.71 

5.38 

1.35 

.693 - 1.972 

A1 2024 

2.785 

5.355 

1.345 

.015 - 1.022 

Ni 

8.86 

4.646 

1.445 

1.009 - 1.491 

Pt 

21.43 

3.646 

1.535 

.322 - 2.718 
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Figure 3.- Wave diagram for material and wave motions occurring along x axis. 
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Figure 4.- Normalized penetration o£ simulating projectile versus fineness 
ratio of corresponding polyethylene cylinder impacting at 11.3 km/sec. 
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Figure 6.- Plaster replica (from ref. 2) of impact crater in 2024-T351 
aluminum target produced by 1/3-fineness-ratio polyethylene cylinder 
at 11.3 km/sec. 
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1 / d =0.151 



v = 7.82 km /sec 


Figure 8.- Impact crater in 2024-T351 aluminum target produced by copper 
simulating projectile at Hugoniot pressure near 1.70 megabars. 
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